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New photometry and astrometry of the isolated neutron 
RX J0720.4-3125 using recent VLT/FORS observations* 
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Since the first optical detection of RX J0720.4-3125 various observations have been performed to determine astrometric 
and photometric data. We present the first detection of the isolated neutron star in the V Bessel filter to study the spectral 
energy distribution and derive a new astrometric position. At ESO Paranal we obtained very deep images with FORS 
1 (three hours exposure time) of RX J0720.4-3125 in V Bessel filter in January 2008. We derive the visual magnitude 
by standard star aperture photometry. Using sophisticated resampling software we correct the images for field distortions. 
Then we derive an updated position and proper motion value by comparing its position with FORS 1 observations of 
December 2000. We calculate a visual magnitude of = 26.81 ± 0.09 mag, which is seven times in excess of what 
is expected from X-ray data, but consistent with the extant U, B and R data. Over about a seven year epoch difference 
we measured a proper motion of /i = 105.1 ± 7.4masyr~^ towards 6 = 296.951° ± 0.0063° (NW) , consistent with 
previous data. 
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1 Introduction 

In deep optical follow-up obs ervations of bright X-ray sour- 
ces from the ROSAT mission dVoges et al seven ther- 
mal emitting isolated neutron stars (INS), exh ibiting no ra- 
dio emission, have been discovered so far (see lHaberlll2QQ4l 
2007 ; van Kerkwiik & Kaplanll2007 for recent reviews). 

After the first optical detection of the brightest one, RX 
J1856.4-3754 (Walter, Wolk, & Neuhiiuser 1996), 



Haberl et all (119971) found the second brightest INS 
RXJ072 0.4-3125. It was detected in the optical B and R 
band bv iKulkarni & van KerkwiikI (Il998h . Using deep ob- 
servations Motch, Zavlin, & Haberl (2003) derived the pro- 
per motion of RX J0720.4-3125 and measured the B mag- 
nitude with VLT/FORS 1. Later Kaplan, van Kerkwijk, 
& Anderson (2007) measured the distance of RX J0720.4- 
3125 ('--^ 360 pc) and updated the proper motion using Hub- 
ble Space Telescope (HST) observations. 

Comparing blackbody fits derived from X-Ray observa- 
tions with the flux m easured in the optical/UV wavelength 
jMotch et al. Il2003h shows that RX J0720.4-3125 exceeds 



its expected optical/UV fl ux by about a n order of magnitude 
jHambaryan et al.ll2009t) . Furthermore iHaberl et all(l2006l) 
detected variations in the X-ray flux. The origin of these 



* Based on observations made with ESO Telescopes at the Paranal Ob- 
servatory under programme ID ()80.D-0135(A). Based on observations 
with XMM-Newton, an ESA Science Mission with instruments and con- 
tributions directly funded by ESA Member states and the USA (NASA). 
** Corresponding author: e-mail: eisen@astro.uni-jena.de 



variations is still controversial (Hohle et al. 2009; Kaplan 
et al. 2007). 

In this paper we measure the V magnitude for the first 
time and update the proper motion and absolute position. 

In Section |2] we discuss the observations and data re- 
duction. Sections|3]and|4]addresses the photometry. In Sec- 
tion |5] we discuss how our newly derived V magnitude fits 
with the investigations of earlier work. Section |6] is dedi- 
cated to our astrometric investigations and in Section [T] we 
summarize our results. 

2 Observations and Data reduction 

In January of 2008 we observed RX J0720.4-3125 for about 
three hours exposure time at ESO Paranal observatory with 
the FOcal Reducer/low dispersion Spectrograph (FORS 1, 
lAppenzeller et alJll998h . installed in Unit Telescope Kue- 
yen of the Very Large Telescope (VLT). Observations were 
carried out in the V Bessel filter. We used a pixel scale of 
0.125"/pixel which leads to a field of view of 4.2' x 4.2'. 13 
single exposures of 900 s were taken over 1.5 nights. 

For astrometric calibration we observed the globular 
cluster 47-Tuc with a pointing at a = 00''22™29.3'* and 
S = — 71'^59'54.3". The observations are summarized in 
Tabled The image of the RXJ0720.4-3 125 field is shown 
in Fig.[T] 

All images were flat fielded and bias subtracted. Since 
FORS 1 got a new detector in 2007 we were restricted to 
absolute world coordinates for astrometric comparison to 
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Table 1 Observation log 



Date 


Instrument 


Filter 


Exp 


Pixel 


Seeing 


Target 




telescope 




s 


(") 


(") 




12 Jan. 2008 


FORSWLT-UTl 


V Bessel 


8.1 • 10^ 


0.125 


1.1 


RXJ0720.4-3125 


13 Jan. 2008 


FORSl/VLT-UTl 


V Bessel 


3.6 • 10^ 


0.125 


0.8 


RXJ0720.4-3125 


12 Jan. 2008 00:42 (UT) 


FORSl/VLT-UTl 


V Bessel 


4 ■ 10 


0.125 


1.6 


47-Tuc 


12 Jan. 2008 00:47 (UT) 


FORSl/VLT-UTl 


V Bessel 


12 


0.125 


1.6 


Landolt field SA 98 


12 Jan. 2008 04:48 (UT) 


FORSWLT-UTl 


V Bessel 


5 


0.125 


1.1 


Landolt field SA 98 




RXJ 0720.4-3125 '.' 



■ ■■■■ 



j 



Fig. 1 Co-added V band images obtained with FORS 1 on Kueyen in January 2008. The circle shows the position of 
RXJ0720.4-3125 as measured in the FORS 1 image of December 2000 dMotch et al. Il2003h . The arrow labeled 01/2008 
points to our detected position of RX J0720.4-3125 . 



earlier observations. Furthermore the new chip is divided 
in two parts which limits the number of reference stars on 
each part. In order to derive an accurate astrometric solution 
we observed the globular cluster 47-Tuc at four dither po- 
sitions. Based on the ESO fits header we derived image, as 
well as world coordinates of e ach object in th e images with 
the Source Extractor (SE, Ber tin & Arnoutsl [l996). These 
object catalogs are fitted with refere nce to the the tw o mi- 
cron all sky survey catalog (2MASS. ICutri et alJ2003h . pro- 
vided by the ViZiR database (Ochsenbein, Bauer, & Mar- 
cout 2000). After the WCS frame is calibrated in that way a 
fifth order polynomial is fitted to the data. This -algorithm 
corrects the rem aining field distortions. We used the soft- 
ware SCAmM ( Benir&2006 ) for these calculations. SCAMP 
stores the calibration information (WCS transformation as 
well as field distortion coefficients) in an external header for 
each image. We extracted the field distortion coefficients. 



Software for Calibrating AstroMetry and Photometi'y 



we derived calibrating the 47-Tuc images and used them 
as input for the calibration procedure of the RXJ0720.4- 
3125 images, which is done by SCAMP as well in the same 
way. Similar, but directly without using a calibration clus- 
ter dieR2£J07204;31 25 images of the year 2000 observed 
by IMotch et al. I (l2003h were calibrated. The field distorion 
map of this calibration is shown as an example in Fig.|2] 

Another software called SWarp ( Bertin et al.ll2002h uses 
the external headers produced by SCAMP to align, resam- 
ple, and co-add the images. A mesh based background sub- 
traction is applied during co-addition and the two parts of 
the chip are stitched together This is illustrated in Fig. [3] 

Typically the semi major axis of the error ellipse of 
2MASS does not exceed 0.3 arcsec. The statistical errors of 
the source detection procedure are much smaller SCAMP 
preselects unambiguous, well exposed sources for the fit- 
ting and treats the uncertainties in a sense. Nevertheless 
it is unnecessary to treat each error source individually. The 
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Instrument Al: distortion map 



I 



Fig. 2 Field distortion map of FORS 1 detector in Decem- 
ber 2000. Color coding goes from red (middle in the image 
and top of the scale) through yellow and green to blue (cor- 
ners of the image and bottom of the scale.) 



resulting measurement plus calibration error is calculated in 
a statistical way as described in Section|6l 



3 Photometry for January 12th, 2008 

To flux-calibrate th e image we ob served the Landolt stan- 
dard star field SA98 (lLandolllll992h in the first night. Since 
the standards were taken only in the first night we used the 
nine images of RX J0720.4-3125 taken in the first night for 
standard star photometry. We co-added them in the same 
way as described above, but only the upper part of the chip 
was used this time to avoid possible systematic photometry 
errors as they may happen while equalizing the background 
of the two parts of the FORS 1 detector In the upper half of 
the chip only the two standard stars 98 556 and 98 557 were 
observed at airmasses 1.656 and 1.112, respectively, (Fig. 

During data analysis it turned out, that there are faint ob- 
jects nearby. Using the Landolt magnitudes for the stars and 
the given errors we first calculated the difference between 
the two magnitudes of the stars as measured by Landolt and 
by us assuming that 



556, Land 



557, Land 



556, inst 



557, inst ; 



where V^„st are our instrumental magnitudes. Our measure- 
ments show that the stars mar ked as c 1 and c2 in Fig. |4] 
were included in the aperture of lLandolti ( Il992h so we took 
them into account accordingly. Aperture phot ometry is done 
using the DAOPHOT package (' StetsonlfigSTi) . 

We measured the instrumental magnitudes of both stan- 
dard stars in each field and calculated the zero point correc- 
tion c and the first order extinction coefficient k. Taking into 
account all sources of error, which is the 3(T error of the in- 
strumental magnitude, the error of the Landolt magnitudes 
and the difference between the values by using the different 



Fig. 3 Illustration how SCAMP co-added the 13 FORS 
1 images (two parts of the chip). We used a dither pat- 
tern of four positions to avoid bad pixels. Diamonds and 
pluses symbolize objects found in each, or at least some of 
the images and used for calculation of the astrometric solu- 
tion while squares symbolize objects not used for calibra- 
tion. The rectangles illustrate the chip shape and the relative 
position. 



Landolt field S98 




Fig. 4 VLT/FORSl image of the Landolt standard stars 
98 556 and 98 557. There are two objects close to the stan- 
dard stars, marked as cl and c2, which contribute flux to 



the standard magnitude measured by iLandolti (Il992h so we 
measured the flux of these objects as well and took their 
contribution into account. 



standard stars we derive 

c = (-21.0832 ±0.0017) mag and 
k = (0.1562 ±0.0058) mag. 

Since these quantities are derived using two different stars, 
their uncertainties include an estimation of the error, which 
is made neglecting the color term. 

The most pronounced error source comes from the faint 
INS magnitude measurement itself. Fig. [T] shows the rim of 
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pixel 

Fig. 5 Projection of the detection of RX J0720.4-3 125 
The neutron star is detected with ^ 6(t. 



an association of bright O-stars to the left of RX J0720.4- 
3125 . In our deep exposures these stars influence the back- 
ground, so we had to choose an aperture smaller than that 
in the standard star image to deal with this problem. The 
INS is detected at ~ Q a. Because of the faintness of the 
source only the top of the PSF has a higher count rate than 
the background hence, for the INS a smaller aperture than 
for the reference stars is possible (Fig. |5]l. We investigated 
the magnitude - aperture dependence of this faint source and 
concluded from our measurements and Fig. |5] that an aper- 
ture radius of seven pixels is sufficient. 

Using the uncertainties of c and k and taking the vari- 
ation of the airmass (Yns — 1-028 ± 0.019) during obser- 
vation into account we finally derive the V magnitude of 
RXJ0720.4-3125tobe 

V = (26.88 ± 0.15) mag, for Jan 12th, 2008. 




Fig. 6 Co-added and background subtracted image used 
for relative photometry. The reference sta r s used are named 
and marked by arrows. iMotch & Haberll (Il998h measured 
their V magnitude with an accuracy of 0.05 mag. The INS 
is marked and the aperture radius (5 pixel, circle) is shown. 



chose smaller inner and outer sky radii compared to the ref- 
erence stars to prevent the neighboring stars from falling 
into the background annulus. 

The difference of the instrumental and the apparent 
magnitudes for each reference star is calculated. The mean 
of this difference gives again the detector zero point, so this 
the magnitude of RX J0720.4-3 125 is determined to be 

V = (26.81 ± 0.09)mag, for Jan 13th, 2008 

which is fully consistent with our other measurement. 

Note, that if the small aperture would have caused any 
loss of flux the magnitude of the INS would be larger in- 
stead. Since the error bar is smaller we use the relative value 
in the following discussion. 



4 Photometry for January 13th, 2008 

Four images were taken in the second night under much bet- 
ter seeing conditions than the first night. Since the most pro- 
nounced source of error in the absolute photometry was the 
measurement of the faint INSs magnitude itself, we decided 
to use these images too, even though no official photomet- 
ric standards were taken that night. We co-added the im- 
ages as before and explained in Section |2] In addition we 
were using the SE to apply a background subtraction. The 
SE uses a mesh based background algorithm with a mesh 
size of 64 pixels. In the resulting image we identified three 
reference st ars near RX J0720.4-3 1 25 with known V-band 
magnitudes (IMotch & Haberllll998l andFig.l6]). 

We used an aperture of 15 pixel radius for the reference 
stars, given by the full width at half maximum (FWHM) 
size of « 4.5 pixel. For the INS we determined the ob- 
ject's size. It turned out that the PSF of the INS vanishes 
below the background at a radius of less than five pixels, so 
we chose the aperture radius accordingly. Furthermore we 



5 The Spectral Energy Distribution of 
RXJ0720.4-3125 

The visual magnitude of RX J0720.4-3125 is consistent 
with an optical excess about an order of magnitude larger 
than the expected flux extrapolated from t he X-r ay spectra 
first reported in Kulkarni & van Kerkwijy (119981) and stud- 
ied in detail by iKaplan et al.l (l2003l) . However, the origin 
of this enigmatic property as well as its relation to the X- 
ray var i ations is sti ll unknown and a m atter of debate, see 
Haberll ( l2007h and iHaberl et al] (|2006|) for a review about 
the INSs and RX J0720.4-3125 in particular. It might be 
possible that the emitting area of the soft X-ray radiation is 
not strictly connected to the source of op tical/UV photons 
dHambarvan et al.ll2009l) . iTriimped (l2005h : Turolla, Zane & 
Drake (2004); Zane, Turolla & Page (2007); Zane, Turolla 

6 Drake (2004) and references therein discuss alternative 
explanations of the optical excess. The radius of the X-Ray 
emitting area is ~ 4.5 km (Haberl et al.. 20061 normalized 
to a distance of 300pc. 
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Table 2 Optical, ground based observations of RX 
shown, in contrast to Fig.|7] where HST observations 
the table. 



J072 0.4-3125. Only meas urements of ground based observations are 
from lKaplan et alJ (l2003h are included. Other references are given in 



Filter 


Mag 




Flux density 
[W/m^/Hz] •lO"-^'' 


•^ccnt 

r ! T 

[A] 


Flux 

[erg s ^ cm ^\ -10 


Refemence 


U 


25.68 ± 


0.17 


9.75till 


3600 




Motch et al. 2003 


B 


26.58 ± 


0.25 


n r aA-1 46 


4300 




Motch et al. 2003 




26.79 ± 


0.20 


r-, r- aA-\ fil 


4300 


r- , -1-1 12 


Motch et al. 2003 




26.44 ± 


0.15 


11.0011:1° 


4300 


7.67ll:J? 


Motch et al. 2003 




26.787 ± 0.040 


'•^"-0.28 


4300 


<: rr+0.21 
J.OO_0.20 


Motch et al. 2003 




26.620 ± 0.050 


q 7q+0.44 
^•^^-0.42 


4300 


6.48l°:lJ 


Motch et al. 2003 


V 


26.81 ± 


0.09 


C £0+0.99 

6.69_o 86 


5500 


3.65l°:r, 


this paper 


R 


26.9 ± 


0.3 


c ,,+1.63 
■'•^^-1.23 


7000 


^■^^-0.53 


Kulkarni & van Kerkwijk 1998 


H 


>23.07 ± 


0.11 


, 7 7^+0.80 
^ '■'J-0.72 


16500 


< 12.29tl:?2 


Posselt et al. 2009 (upper limit) 



For comparison to other optical/UV magnitudes we 
show in Fig. |7]the flux expected from a black body model 
(XMM-Newton EPIC-pn). The effective temperature of 
RX J0720.4 - 3 1 25 is ch anging on time scal es of years, see 
Hohle et al.l (l2009l) and lHaberl etaP (l2006h . 

We plot in Fig. |7] the spectra obtained with 
XMM-Newton EPIC-pn with lowest (orbit 0078) and high- 
est (orbit 0815) temperatures (86.5 ± 0.4 eV and 94.6 ± 
0.5 eV, respectively) measured yet. The values were obtain- 
ed by fitting the spectra with standard software xspecl2 
(Dorman, Arnaud, & Gordon 2003) in the energy range of 
0.16-1.5 keV by applying a black body model with addi- 
tive gaussian absorption line (line energy = 301 ± 3eV, 
cr = 77 ± 2 eV) wifli absor ption due to die IS M {Nh = 
1.04±0.02-1020 cm-2), see 



Hohle et alj ( l2009t) for details. 



6 Astrometry 



For the longe st possible epoch di fference we used the B 



band image of iMotch et al. I (1200 3h from the year 2000 and 
compared the position of all objects with the positions in 
our image from 2008. The objects in the two images are as- 
signed to each other and the change of position is calculated 
(Fig.®. 

Afterwards the change of position undergoes a Kolmo- 
gorov-Smimov 2-sample test comparing the absolute posi- 
tion changes of the objects to a simulated Rayleigh distri- 
bution. The test is performed comparing the cumulative dis- 
tribution functions (CDF) of the test sample and the data. 
If the test fails objects lying outside la from the mean of 
the distribution are excluded. The mean and standard devia- 
tion of the distribution are recalculated and the motion of all 
objects is shifted by that mean to reduce the last remaining 
systematical errors. This is repeated until the test succeeds. 
What is left is a sample of Rayleigh distributed background 




frequency [Hz] 



Fig. 7 Optical/UV flux from R to UV magnitudes (our 
new V-band data point is marked as circle) compared to X- 
Ray spectra obtained from XMM-Newton EPIC-pn (right) 
performed in full frame mode with thin filter The two spec- 
tra of RXJ0720.4-3125 are taken from orbit 0078 (grey 
dots) and 0815 (black dots) with effective temperatures of 
86.5±0.4eV and 94.6±0.5eV, respectively. The solid fine 
marks the best fit model for rev. 0078 (including interstel- 
lar extinction), while the dashed lines corresponding to the 
black bodies for both spectra without any absorption. The 
broad absorption feature near 300eV (=7.25 • 10^*^ Hz) can 
be seen in the spectrum of revolution 0815 (black dots). 
The Optical magnitudes are sumarized in table |2l the flux 
in the H band is an upper limit. 



stars, without any systematic effects. The standard devia- 
tion (TBack of these stars could be seen as total calibration er- 
ror Since the NS is much fainter than the background stars, 
the positional error of its photo-center is relatively large. 
We added the intrinsic proper motion error A/iintr (which 
is again calculated by the position error in both epochs and 
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Our proper motion (/i = 105.1ib7.4 mas) co nfirms the mea 



surement obtained by iKaplan et al. I (l2007h with the HST 



10 



12 



14 



Fig. 9 Cumulative distribution function (CDF) of the 
Rayleigh distributed background stars in the RXJ0720.4- 
3125 field, compared with the CDF of a synthetic Rayleigh 
distribution. The comparison shows good agreement, prov- 
ing that there are no systematical errors left in the calcula- 
tion. 

the epoch difference) of the NS in both epochs to that error, 
representing the full error of the neutron stars motion as 



A/iNS = \J (A^intr)^ + (o-back)^- 

The mean intrinsic error is already included in tiBack but 
is larger for the barely visible NS than for the background 
objects. 

We derive a proper motion value 
(-93.2 ± 5.4) mas/yr 
(48.6 ± 5.1) mas/yr 
(105.1 ± 7.4) mas/yr. 
B = (296.951 ±0.0063)° 
This proper mot ion is co ns istent wi th previously pu blishe d 



(107.8 ±1.2 mas) with a difference of less than 3 mas/yr. 

Our V magnitude (V = 26.81 ± 0.09 ma g) is compat- 
ible with previous observations from Motch et al. I (2003) 
(B = 26.620 ±0.050 mag) and from (Kulkarni & van Kerk- 
wijk (1998)) (R = 26.9 ± 0.3 mag). These data can be fit- 
ted by the emergent spectrum of the photoi onized plasma, 
surrounding the NS (iHambarvan et al 1 l2009h . This scenario 
would also explain the power-l aw component of the "opti- 
cal excess" (iKaplan et al.ll2003h . which is confirmed by our 
data. 

For further investigations correlated observations in X- 
ray and the optical/UV are needed to collect enough data to 
investigate the nature of the flux variations of RX J0720.4- 
3125 and to understand this isolated neutron star. 

We also note, that there are no additional faint, fast mov- 
ing objects detected in this field (by comparing the year 
2000 observations and our image), hence no additional sim- 
ilar NS in this field. 
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values (Kapla n etal.1 120071: iMotch et al. I l2003b . References 



RX J0720.4-3 125 changed its position from December 2000 
( MJD = 51902 days) to January 2008 ( MJD = 54477 
days) from 



MJD = 51902 days 
1^ 20™ 24.976" 
± 14.15 mas 
MJD = 54477 days 
1^ 20" 24.926" 
± 19.78 mas 



-SI'' 25' 50.105" 
± 12.82 mas 

-31'^ 25' 49.776" 
± 19.30 mas, 



to 



Since the total exposure time in the year 2000 image was 
longer and the observation wavelength was shorter, our 
newly reduced position measurement of the year 2000 is 
more precise and the signal to noise ratio is better 

7 Summary 

Taking a deep optical image of RX J0720.4-3125 we were 
able to determine the V magnitude and updated the position. 
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Fig. 8 Detected objects within a radius of 1' around the position of RX J0720.4-3125 . The lines indicate the change 
of position between the two images of MJD 51902 and MJD 54477 of the objects, scaled up by a factor of 20. While 
RX J0720.4-3125 (middle) is obviously moving, most of the other objects are almost standing still, demonstrating the high 
accuracy of the calibration. Other "fast moving" objects are false detections or wrong assignments of objects between the 
two images. 
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